The Astrophysical Journal, 679, 2008 May 20, in press 

Preprint typeset using I4T]^X style emulateapj v. 03/07/07 



00 

O 
O 

(N 



6 



> 



O 

oo 
O 



IMPROVED CONSTRAINTS ON THE PREFERENTIAL HEATING AND ACCELERATION OF OXYGEN IONS IN THE 

EXTENDED SOLAR CORONA 

Steven R. Cranmer, Alexander V. Panasyuk, and John L. Kohl 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
Submitted 2007 November 27; accepted 2008 January 29 
Draft version February 4, 2008 

ABSTRACT 

We present a detailed analysis of oxygen ion velocity distributions in the extended solar corona, based on 
observations made with the Ultraviolet Coronagraph Spectrometer (UVCS) on the SOHO spacecraft. Polar 
coronal holes at solar minimum are known to exhibit broad line widths and unusual intensity ratios of the O VI 
AA1032, 1037 emission line doublet. The traditional interpretation of these features has been that oxygen ions 
have a strong temperature anisotropy, with the temperature perpendicular to the magnetic field being much 
larger than the temperature parallel to the field. However, recent work by Raouafi and Solanki suggested that 
it may be possible to model the observations using an isotropic velocity distribution. In this paper we analyze 
an expanded data set to show that the original interpretation of an anisotropic distribution is the only one that 
is fully consistent with the observations. It is necessary to search the full range of ion plasma parameters to 
determine the values with the highest probability of agreement with the UVCS data. The derived ion outflow 
speeds and perpendicular kinetic temperatures are consistent with earlier results, and there continues to be 
strong evidence for preferential ion heating and acceleration with respect to hydrogen. At heliocentric heights 
above 2. 1 solar radii, every UVCS data point is more consistent with an anisotropic distribution than with an 
isotropic distribution. At heights above 3 solar radii, the exact probability of isotropy depends on the electron 
density chosen to simulate the line-of-sight distribution of O VI emissivity. The most realistic electron densities 
(which decrease steeply from 3 to 6 solar radii) produce the lowest probabilities of isotropy and most-probable 
temperature anisotropy ratios that exceed 10. We also use UVCS O VI absolute intensities to compute the 
frozen-in O s+ ion concentration in the extended corona; the resulting range of values is roughly consistent with 
recent downward revisions in the oxygen abundance. 

Subject headings: line: profiles — plasmas — solar wind — Sun: corona — Sun: UV radiation — techniques: 
spectroscopic 



1. INTRODUCTION 

The physical processes that heat the solar corona and ac- 
celerate the solar wind are not yet understood completely. 
In order to construct and test theoretical models, there must 
exist accurate measurements of relevant plasma parameters 
in the regions that are being heated and accelerated. In 
the low-density, open-field regions that reach into interplan- 
etary space, the number of plasma parameters that need to 
be measured increases because the plasma begins to become 
collisionless and individual particle species (e.g., protons, 
electrons, and heavy ions) can exhibit different properties. 
Such differences in particle velocity distributions are valu- 
able probes of "microscopic" processes of heating and accel- 
eration. The Ultraviolet Coronagraph Spectrometer (UVCS) 
operating aboard the Solar and Heliospheric Observatory 
(SOHO) spacecraft has measured these properties for a va- 
riety of open-field regions in the extended corona (Kohl et al. 
1995, 1997, 2006). 

In this paper we focus on UVCS observations of heavy 
ion emission lines (specifically O VI AA1032, 1037) in po- 
lar coronal holes at solar minimum. One main goal is to re- 
solve a recent question that has arisen regarding the existence 
of anisotropic ion temperatures in polar coronal holes. Sev- 
eral prior analyses of UVCS data have concluded that there 
must be both intense preferential heating of the O s+ ions, in 
comparison to hydrogen, and a strong field-aligned anisotropy 
with a much larger temperature in the direction perpendicular 
to the magnetic field than in the parallel direction (see, e.g., 
Kohl et al. 1997, 1998; Li et al. 1998; Cranmer et al. 1999; 



Antonucci et al. 2000; Zangrilli et al. 2002; Antonucci 2006; 
Telloni et al. 2007). However, Raouafi & Solanki (2004, 
2006) and Raouafi et al. (2007) have reported that there may 
not be a compelling need for O s+ anisotropy depending on the 
assumptions made about the other plasma properties of the 
coronal hole (e.g., electron density). 

The determination of O s+ preferential heating, preferential 
acceleration, and temperature anisotropy has spurred a great 
deal of theoretical work (see reviews by Hollweg & Isenberg 
2002; Cranmer 2002a; Marsch 2005; Kohl et al. 2006). It 
is thus important to resolve the question of whether these 
plasma properties are definitively present on the basis of the 
UVCS/SOHO observations. In this paper, we attempt to ana- 
lyze all possible combinations of O s+ properties (number den- 
sity, outflow speed, parallel temperature, and perpendicular 
temperature) with the full effects of the extended line of sight 
(LOS) taken into account. The applicability of any particu- 
lar combination of ion properties is evaluated by computing 
a quantitative probability of agreement between the modeled 
set of emission lines and a given observation. Preliminary re- 
sults from this work were presented by Cranmer et al. (2005) 
and Kohl et al. (2006). 

The original UVCS results of preferential ion heating and 
acceleration — as well as strong ion temperature anisotropy 
(T± 3> 7j|) — were somewhat surprising, but these extreme de- 
partures from thermal equilibrium are qualitatively similar to 
conditions that have been measured for decades in high-speed 
streams in the heliosphere. At their closest approaches to 
the Sun (~ 0.3 AU), the Helios probes measured substantial 
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proton temperature anisotropics with 7j_ > 7]| (Marsch et al. 
1982; Feldman & Marsch 1997). In the fast wind, most ion 
species also appear to flow faster than the protons by about 
an Alfven speed (V A ), and this velocity difference decreases 
with increasing radius and decreasing proton flow velocity 
(e.g., Hefti et al. 1998; Reisenfeld et al. 2001). The tem- 
peratures of heavy ions are significantly larger than proton 
and electron core temperatures. In the highest-speed wind 
streams, ion temperatures exceed simple mass proportional- 
ity with protons (i.e., heavier ions have larger most-probable 
speeds), with (T ion /T p ) > (m ion /m p ), for m ioB > m p (e.g., Col- 
lier et al. 1996). UVCS provided the first evidence that these 
plasma properties are already present near the Sun. 

The outline of this paper is as follows. In § 2 we present an 
expanded collection of WCS/SOHO observational data that 
is used to determine the 5+ ion properties. § 3 outlines the 
procedure we have developed to produce empirical models 
of the plasma conditions in polar coronal holes and to com- 
pute the probability of agreement between any given set of 
ion properties and the observations. The resulting ranges of 
ion properties that are consistent with the UVCS observations 
are presented in § 4 along with, in our view, a resolution of 
the controversy regarding the oxygen temperature anisotropy. 
Finally, § 5 gives a summary of the major results of this paper 
and a discussion of the implications these results may have on 
theoretical models of coronal heating and solar wind acceler- 
ation. 

2. OBSERVATIONS 

The UVCS instrument contains three reflecting telescopes 
that feed two ultraviolet toric-grating spectrometers and one 
visible light polarimeter (Kohl et al. 1995, 1997). Light from 
the bright solar disk is blocked by external and internal occul- 
ters that have the same linear geometry as the spectrometer 
slits. The slits are oriented in the direction tangent to the solar 
limb. They can be positioned in heliocentric radius r any- 
where between about 1.4 and 10 solar radii (Rq) and rotated 
around the Sun in position angle. The slit length projected 
on the sky is 40', or approximately 2.5 R Q in the corona, and 
the slit width can be adjusted to optimize the desired spectral 
resolution and count rate. 

The UVCS data discussed in this paper consist of a large 
ensemble of observations of polar coronal holes from the last 
solar minimum (1996-1997). The solar magnetic field is ob- 
served to exist in a nearly axisymmetric configuration at solar 
minimum, with open field lines emerging from the north and 
south polar regions and expanding superradially to fill a large 
fraction of the heliospheric volume. The plasma properties 
in polar coronal holes remain reasonably constant in the year 
or two around solar minimum (see, e.g., Kohl et al. 2006), 
so we assemble the data over this time into a single function 
of radius. In this paper, we limit ourselves to the analysis of 
observations of the O VI AA1032, 1037 emission line doublet 
in these polar regions. The relevant UVCS observations are 
taken from the following three sources. 

1. The empirical model study of Kohl et al. (1998) and 
Cranmer et al. (1999) covered the period between 1996 
November and 1997 April and took the north and south 
polar coronal hole properties to be similar enough to 
treat them together. 

2. A detailed analysis of the north polar coronal hole 
by Antonucci et al. (2000) coincided with the second 



SOHO joint observing program (JOP 2) on 1996 May 
21. 

3. We searched the WCS/SOHO archive for any other 
north or south polar hole observations having sufficient 
count-rate statistics to be able to measure the O VI line 
widths at radii above 2 Rq. A total of 14 new or rean- 
alyzed data points were identified between 1996 June 
and 1997 July. 

The remainder of this section describes the data reduction for 
the third group of new data points. Table 1 provides details of 
these 14 measurements. 

The criteria for identifying new UVCS data were as fol- 
lows. We adopted a time period from 1996 April, the be- 
ginning of primary science operations, to 1998 January, after 
which the new cycle's activity began to rise and high-latitude 
streamers appeared regularly to signal the end of true solar 
minimum. Only measurements of the O VI lines above the 
poles (i.e., position angles within ±15° of the north or south 
poles) and at heights above 2 Rq were sought. 1 Prior expe- 
rience with the count rates at large heights in coronal holes 
refined the search further to use measurements only with rel- 
atively long exposure times (see Table 1) to gather sufficient 
statistics to measure the line widths. There were two observa- 
tions that appeared initially to satisfy the above criteria (1997 
April 15, at 4.14 R Q , and 1997 July 2, at 3.10 R Q ), but they 
were not used because the count rate statistics were inade- 
quate for reliable line widths. The only point of overlap be- 
tween the data in Table 1 and prior analyses (e.g., Cranmer 
et al. 1999) concerns the end of the month-long study of the 
north polar coronal hole in 1997 January (at r k, 3Rq). These 
data were reanalyzed with a different line fitting technique 
and an improved UVCS pointing correction; the computed 
line widths are similar to those presented by Cranmer et al. 
(1999) and the intensities are given here for the first time. For 
completeness, though, both the old and new data points are 
kept in the full ensemble of O VI data used below. 

To achieve the lowest uncertainties in the determinations of 
the O VI AA1032, 1037 intensities and line widths, we typ- 
ically integrated over 15' to 30' along the slit (see Table 1). 
This corresponds to ±0.5-1 Rq on either side of the north- 
south axis. The use of such large areas implies that narrow 
flux-tube structures such as dense polar plumes and the less- 
dense interplume regions were not resolved. As long as all 
steps of the analysis remain consistent with such a coarsely 
averaged state (e.g., the use of a similarly averaged electron 
density), though, this need not be a problem. The derived 
plasma properties thus describe the average conditions inside 
coronal holes at solar minimum and do not address differences 
between plumes and interplume regions. 

Details concerning the analysis of UVCS data are given by 
Gardner et al. (1996, 2000, 2002) and Kohl et al. (1997, 1999, 
2006). The UVCS Data Analysis Software (DAS) was used 
to remove image distortion and to calibrate the data in wave- 
length and intensity. The coronal line profiles are broadened 
by various instrumental effects. The optical point spread func- 
tion of the spectrometer depends on the slit width used (with 
270.3 fim corresponding to 1 A in the spectrum), the on-board 

1 Below r 2Rq, the existing data appear to be adequate, uncertainties 
are low, and there is not much of an intrinsic spread in the intensities and 
line widths as a function of height. Also, earlier analyses did not show that 
collisionless effects (ion temperature anisotropics, preferential ion heating, or 
differential flow) became strong until above this height. 
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TABLE 1 

Newly Analyzed Polar Coronal Hole Data: 1996- 
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Mart Date, U 1 lime 


Obs. Height 


Mit Length 


Mit Width 


Exposure 


Vl/e 
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itot, 1U32 A line 






(arcmin) 
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Time (hr) 


(kms 4 ) 




(10 6 phots 1 cm 2 sr ') 


1 QQ£ Tun 91 1 ft-^s 
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9 A7 
Z.U / 


1 J.o ) 


/ 9 


Q 1 
y. 1 


TAT -4-917 
JO J IE Zl . / 


9 no -t- n A 


JZ.Z IE O.o 


LyytJ jcp Zy, 19. ZD 


9 fi8 
Z.Uo 


1 / . J ^1N ) 


7S 

/ 9 


9.D 


AIT -4- 1 9 A 


1 Q7 _t_ ft 1 A 
1 .0 / ZH U. 1^+ 


/I a 9 _l_ 7 
J.Z xo./ 


1996 Nov 07, 23:38 


2.07 


29.7 (N) 


75 


7.6 


409 ± 18.8 


1.78±0.2 


44.6±9.1 


1996 Nov 10, 06:30 


3.00 


25.7 (N) 


350 


9.4 


475 ± 27.0 


1.16±0.18 


2.08 ±0.44 


1996 Nov 10, 15:55 


2.56 


25.9 (N) 


350 


10.8 


505 ± 30.5 


1.17±0.17 


4.83 ±1.0 


1996 Nov 16, 16:56 


2.17 


29.6 (N) 


340 


9.5 


417 ±7.43 


1.485 ±0.06 


25.1 ±5.0 


1997 Jan 05,21:00 


3.07 


28.0 (N) 


100 


17.6 


690 ±87.2 


0.957 ±0.45 


2.02 ±0.63 


1997 Jan 10, 14:54 


3.08 


20.8 (N) 


150 


68.8 


686 ±38.7 


1.23 ±0.3 


3.37 ±0.76 


1997 Jan 24, 16:03 


3.08 


20.6 (N) 


150 


70.8 


594 ±41.6 


1.01 ±0.3 


1.66 ±0.40 


1997 Mar 09, 18:00 


2.57 


19.0 (S) 


300 


8.9 


500 ± 22.0 


1.15 ±0.12 


5.66±1.2 


1997 Jun 04, 16:49 


2.56 


18.8 (N) 


300 


9.0 


527 ±33.9 


0.958 ±0.15 


4. 14 ±0.90 


1997 Jun 08, 20:15 


3.10 


18.7 (N) 


300 


8.3 


645 ±65.6 


1.33 ±0.44 


1.33 ±0.33 


1997 Jul 01, 15:50 


2.56 


17.2 (N) 


342 


9.8 


534 ±56.8 


0.988 ±0.26 


3.26 ±0.80 


1997 Jul 04, 16:45 


3.63 


18.9 (N) 


342 


29.5 


451 ±47.9 


1.49 ±0.43 


0.559 ±0.13 



' Data were integrated over the specified slit length. The slit was oriented tangent to either the north (N) or the south (S) heliographic pole, as indicated. 



data binning, the exposed mirror area, and the intrinsic quan- 
tization error of the detector. This broadening is taken into 
account by adjusting the line widths of Gaussian fits to the 
coronal components of the data; the data points themselves 
are not corrected. Tests have shown that the coronal line width 
can be recovered accurately even when the total instrumental 
width is within about a factor of two of the width of the coro- 
nal component. Instrument-scattered stray light from the solar 
disk is modeled as an additional narrow Gaussian component 
with an intensity and profile shape constrained by the known 
stray light properties of the instrument. 

The analysis of the O VI emission line doublet involves four 
basic observable quantities: the total intensities of the two 
lines and their 1/e Gaussian half-widths AAi/ e . The latter 
quantities are typically expressed in Doppler velocity units as 
V\/ e = cAA^g/Ao, where Ao is the rest wavelength of the line 
and c is the speed of light. Rather than give the two total in- 
tensities, Table 1 provides the total intensity I tot of the O VI 
A 1032 line and the ratio K of the A 1032 to the A 1037 inten- 
sities. The uncertainties given in Table 1 take account of both 
Poisson count-rate statistics and the fact that the various in- 
strumental corrections are known only to finite levels of pre- 
cision. Note that the ratio 1Z does not depend on the absolute 
intensity calibration of the instrument. 

UVCS/SOHO has not been able to resolve any departures 
from Gaussian shapes for the O VI lines in large polar coronal 
holes, so the profiles are described by just the one parameter 
V\/ e . For the measurements given in Table 1, we performed 
the line fitting by constraining the coronal components of the 
A 1032 and A 1037 lines to have the same width. Thus, the 
V\/ e values given in Table 1 are formally a weighted mean be- 
tween the two components. This is done mainly to lower the 
statistical uncertainties but there is some observational justi- 
fication for assuming that the two components have the same 
width. In situations where the count rates are high, it is diffi- 
cult to see any significant or systematic difference between the 
line widths of the two components. There are various reasons 
why they may be different from one another in some regions 
(e.g., Cranmer 2001 ; Morgan & Habbal 2004), but more work 
needs to be done to identify such subtle effects. 

Figure 1 displays the combined ensemble of old and new 
UVCS O VI data for the three main observables: the line 
width V\j e , the dimensionless intensity ratio 1Z, and the abso- 
lute (line-integrated) intensity of the O VI A 1032 line. There 
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FIG. 1. — Collected UVCS polar coronal hole measurements of (a) O VI 
line widths V l/e , (b) ratio of O VI A1032 to O VI A1037 intensities, and (c) 
O VI A1032 line-integrated intensities, with symbols specifying the sources 
of the data (see labels for references). Error bars denote ±1<t observational 
uncertainties. Also shown (dotted lines) are the parameterized fits given by 
Cranmer et al. (1999). SEE LAST PAGE OF PAPER FOR LARGER 
VERSION. 

are a total of 53 separate data points from the three sources 
discussed above, but not all of these points have all three of 
the main quantities: there are 50 values of V\u, 52 values of 1Z 
(with only 49 cases where both V\/ e and 1Z exist for the same 
measurement), and 44 values of / tot . This relative paucity of 
data illustrates the difficulties of measuring the plasma param- 
eters at large heights in polar coronal holes. 

In general, the radial dependences of the O VI quantities in 
Figure 1 are similar to those given by Kohl et al. (1998), Cran- 
mer et al. (1999), and Antonucci et al. (2000). There exists a 
reasonably large spread in the V\/ e values in Figure la above 
r s=a 2.5Rq. This spread exceeds the magnitude of the ±lcr 
uncertainty limits for the individual measurements, and thus 
seems to indicate that there is an intrinsic variability (possi- 
bly temporal) of the O s+ plasma conditions in polar coronal 
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holes above heights where the ions become collisionless. It 
is possible that polar plumes and interplume regions become 
collisionless over different ranges of radius, and thus prefer- 
ential ion heating mechanisms may begin to broaden the O VI 
lines at different rates in the two regions. The observed varia- 
tion in line width may thus depend on the relative concentra- 
tions of plume and interplume plasma along the line of sight 
at different observation times. 

3. EMPIRICAL MODEL PROCEDURE 

The observable properties of the O VI line doublet depend 
on a nontrivial combination of various O s+ plasma parame- 
ters, as well as electron parameters, integrated along the opti- 
cally thin line of sight. In general, then, it is not possible to de- 
rive accurate and self-consistent plasma parameters via a sim- 
ple "inversion" from the line widths and intensities. Rather, 
one must build up a so-called empirical model of the coro- 
nal hole — with the 5+ velocity distribution and other prop- 
erties as free parameters — and synthesize trial line profiles. 
After some procedure of varying the coronal parameters to 
achieve agreement between the synthesized line profiles and 
the observations, the self-consistent empirical model of the 
ion properties can be considered complete. This technique is 
closely related to forward modeling approaches being used in 
other areas of solar physics (e.g., Judge & Mcintosh 1999). 

The use of the term "empirical model" has resulted in a 
bit of confusion regarding what assumptions are embedded 
in the derived plasma parameters. We emphasize that the 
empirical models described here do not specify the physi- 
cal processes that maintain the coronal plasma in its assumed 
steady state. Thus, there is no explicit dependence on "the- 
oretical" concepts such as coronal heating and acceleration 
mechanisms, waves and turbulent motions, or magnetohydro- 
dynamics (MHD), within the empirical models. The derived 
5+ plasma parameters depend on only the observations and 
on well-established theory such as the radiative transfer inher- 
ent in the line-formation process. 

In this section we summarize the forward modeling of O VI 
line profiles for an arbitrary set of coronal parameters (§ 3.1), 
then we describe how these parameters are specified and var- 
ied to produce various empirical model grids (§ 3.2). Finally, 
we present a new method of computing the probability of 
agreement between a given empirical model and the obser- 
vations (§ 3.3), such that no regions of the possible solution 
space are neglected. 

3.1. Forward Modeling 

The O VI line emission in coronal holes comes from two 
sources of comparable magnitude: (1) collisional electron im- 
pact excitation followed by radiative decay, and (2) resonant 
scattering of photons that originate on the bright solar disk. 
The emergent specific intensity of an emission line from an 
optically thin corona is given by 

/+oo 
dx(ff + j™) , (1) 
oo 

where x is the coordinate direction along the observer's line 
of sight (LOS) and j c ° n and f£ s are the collisionally excited 
and resonantly scattered line emissivities, respectively. We 
neglect the relatively weak UV continuum and the Thomson 
electron-scattered components of the spectral lines in ques- 
tion. At a given point in the three-dimensional coronal hole 
volume, the line emissivities are specified by 

ffj^ qn{ j e) n e n^ v (2 ) 



j^ = ^B 12 n 1 J °d V '^R(u',n',v,n)IA&') (3) 

(see, e.g., Mihalas 1978). Here, vq is the rest-frame line- 
center frequency, qu is the collision rate per particle for the 
transition between atomic levels 1 and 2, n\ is the number 
density in the lower level of the atom or ion of interest (here, 
the 2s 2 Sx/2 state of O s+ ), and B12 is the Einstein absorption 
rate of the transition. The emission profile <\> v is assumed to 
be Gaussian. The scattering redistribution function R takes 
the incident frequency v' and photon direction vector n' and 
transforms it into the observed frequency v along the LOS 
direction n. 

The profile <p u and the redistribution function R contain 
the main dependences on the properties of the ion velocity 
distribution. We allow for the possibility of an anisotropic 
5+ velocity distribution by using a bi-Maxwellian function 
(e.g., Whang 1971), with the parallel and perpendicular axes 
oriented arbitrarily with respect to the radial direction in the 
corona; see § 3.2. The emissivity profiles along the LOS are 
modeled with the full effects of the bi-Maxwellian velocity 
distribution and the projected components of the bulk outflow 
speed along the n' and n directions. For the polar coronal hole 
measurements being modeled here, we define a Cartesian co- 
ordinate system for which the LOS direction is denoted x and 
the north-south polar axis of the Sun is z. The other coordi- 
nate y is set to zero. General expressions for the emissivities 
are given in various levels of detail by Withbroe et al. (1982), 
Noci et al. (1987), Allen et al. (1998), Cranmer (1998), Li et 
al. (1998), Noci & Maccari (1999), Kohl et al. (2006), and 
Akinari (2007). 

The resonantly scattered components depend sensitively on 
the intensity profiles incident from the solar disk (/„')■ As in 
Cranmer et al. (1999), we used empirically derived Gaussian 
profiles with total intensities measured on the disk by UVCS 
at solar minimum (Raymond et al. 1997). The adopted O VI 
A1032 (1031.93 A) disk intensity is 1.94 x 10 13 photons s _1 
cm" 2 sr" 1 , and the total intensities of the O VI A1037 (1037.62 
A), C II A1037.02, and C II A1036.34 disk lines are 0.500, 
0.214, and 0.171 times the A1032 intensity, respectively. We 
used the profile widths as given by Noci et al. (1987); see also 
the comparative tables of Gabriel et al. (2003) and Raouafi & 
Solanki (2004). 

The collisional components depend on how the collision 
rate q\2 varies with electron temperature T e . We kept the 
same tabulated values as were used by Raymond et al. (1997) 
and Cranmer et al. (1999). For completeness, we give a fit to 
qn(T e ) for the O VI A1032 transition: 

loggia) = -0.221 17r 2 + 2.4565r- 14.695 , (4) 

where t = log 10 T e , and T e and qi2 are given in units of K and 
cm 3 s" 1 respectively. This expression is valid to within about 
± 2% over the range 5.3 < t < 6.3. The collision rate for the 
O VI A 1037 line is half of that of the O VI A 1032 line. 

The numerical code that synthesizes line profiles by nu- 
merically integrating equations (HJ-© is essentially the same 
as the one used by Cranmer et al. (1999). The integrations 
over x and v' have been simplified by replacing the adaptive 
Romberg method by fixed grids, with spacings that have been 
adjusted to minimize both numerical discretization errors and 
run time. The LOS integration was performed in steps of 0.1 
Rq from -15 to +15 Rq along the x axis. The incident fre- 
quency grid corresponds to a wavelength grid with a spacing 
of 0.03 A in A'. These step sizes were verified to give accurate 
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results by halving the step sizes and obtaining the same results 
to within a desired precision. We integrated over the solid an- 
gle of the solar disk (dQ! = sin 6' d9' d(j>') by Gauss-Legendre 
quadrature in 0' and equally spaced trapezoidal quadrature in 
(/)'. The solar disk was assumed to be uniformly bright. 

3.2. Parameter Selection for Line Synthesis 

For the O VI doublet, there are three primary observables 
(/tot of A1032, V\/ e , and 1Z) that depend on the LOS distribu- 
tions of four "unknown" quantities as well as a longer list of 
quantities that can be considered to be known independently 
of the UVCS observations. The four unknowns are the ion 
fraction (essentially n\/n e ), the O s+ bulk outflow speed along 
the magnetic field (M;ii), and the parallel and perpendicular 
O s+ kinetic temperatures (Tju and Tj±). The known quanti- 
ties include the electron density n e , the electron temperature 
T e , the incident intensity from the solar disk, and the overall 
magnetic geometry of the coronal hole (i.e., how to compute 
"parallel" and "perpendicular" at any point along the LOS). 
Note that both emissivities (eqs. ED-G)) depend linearly on 
the ion fraction, so that the total intensity 7 tot can be used as 
a straightforward diagnostic of this quantity after the other 
parameters have been determined. The line widths and in- 
tensity ratios do not depend on the ion fraction. This leaves 
two observables (V] i e and 1Z) to specify the values of three ion 
quantities (Ujii, TJii, and Tj±), Although this system is formally 
underdetermined, we can nonetheless put some firm limits on 
the ranges of these quantities and compute the most probable 
values. 

Below, the three 5+ velocity distribution parameters are 
discussed in § 3.2.1 and the other "known" parameters are 
discussed in § 3.2.2. 

3.2.1. Ionized Oxygen Parameters 

We treat the three unknown ion quantities as free parame- 
ters that are varied independently of one another. Other em- 
pirical modeling efforts (e.g., Cranmer et al. 1999; Antonucci 
et al. 2000; Raouafi & Solanki 2004, 2006) have tended to 
use some form of iterative refinement; i.e., they started with 
a specific set of initial conditions and assumptions, and they 
varied some parameters — and kept others fixed — to find the 
most probable values of Kjii, Tj-n, and T,±. The initial esti- 
mates tended to utilize the fact that the line widths are most 
sensitive to Tj±, whereas the line ratios depend mainly on the 
effect of Doppler dimming (and Doppler pumping from the 
C II solar disk lines) and thus are sensitive mainly to the par- 
allel velocity distribution (Ufit and 7jii). These iterative proce- 
dures contain the inherent possibility that some regions of the 
parameter space could be neglected, and thus possibly valid 
solutions could be ignored. In this paper we search the entire 
parameter space by constructing a three-dimensional "data 
cube" which contains all possible combinations of the three 
parameters. 

The three axes of the data cube were chosen to be 2^11, 7}j_, 
and the anisotropy ratio T,± /Tq. The modeled ranges of these 
quantities were made as wide as possible in order to avoid 
missing possibly relevant regions of parameter space. The 
outflow speed um was varied between and 1000 km s" 1 us- 
ing a linearly spaced grid. The perpendicular kinetic tem- 
perature 7}j_ was varied logarithmically between 5 x 10 5 and 
10 9 K. The anisotropy ratio was varied logarithmically be- 
tween 0.1 and 100. There were 50 values of each parameter 
along the three axes of the data cube, and we synthesized 12 



wavelengths — spaced linearly between the line center and 2.7 
A redward of line center — for both O VI lines. Thus, a data 
cube constructed for a specific height in the corona (z) con- 
sisted of 3 x 10 6 (50 3 x 24) individual LOS integrations. 

For each point in a data cube, the scalar values of unt and 
7}j_ were assumed to be those in the plane of the sky (i.e., 
x = 0). For other points along the LOS, the models used 
slightly larger values that are consistent with an assumed ra- 
dial increase in both parameters. Mass flux conservation — 
using the modeled n e (r) and flux tube geometry — was used to 
specify the radial increase in ua\ along the LOS. Earlier em- 
pirical modeling results (specifically, eq. [28] of Cranmer et 
al. 1999) were used to specify the radial increase in T,±. The 
modeled anisotropy ratio 7]j_ /7}ii was assumed to remain con- 
stant along the LOS. It is important to note that the modeled 
radial increases in Uj\\ and T,± were always taken to be relative 
to the plane-of-sky values that were varied freely throughout 
each data cube. Thus, there is no a priori reason for the re- 
sulting most-probable values of these parameters (determined 
via comparisons with observations over a range of heights z) 
to exhibit similar radial increases. 2 

We note that the kinetic temperature quantities Tj± and 7}|| 
may describe some combination of "thermal" microscopic 
motions and any unresolved bulk motions due to waves or tur- 
bulence. Thus, there is a further step of interpretation required 
after the most likely values of these kinetic temperatures have 
been derived from the empirical modeling process. Making 
a definitive separation between the thermal and nonthermal 
components of these temperatures is beyond the scope of this 
paper. However, we can make some qualitative comments on 
the likely ranges of magnitude of these two components based 
on recent theoretical models of Alfven waves in coronal holes; 
see §§4.2 and 4.3. 

Finally, the 5+ ion fraction n\/n e was kept at a constant 
(and arbitrary) value in all of the models. Comparisons be- 
tween the observed and synthesized total intensities were used 
to derive measurements of this ion fraction in the polar coro- 
nal holes; see § 4.5. 

3.2.2. Electron and Flux Tube Parameters 

The three main "known" parameters that are explored in the 
models shown below (but kept constant over each data cube) 
are the electron density n e (r), electron temperature T e (r), and 
the macroscopic flux-tube geometry of the coronal hole. Any 
other parameters that could be varied — e.g., the disk intensi- 
ties of the O VI and C II lines — were kept fixed at the values 
given above in § 3.1. 

Because one main purpose of this paper is to determine why 
the results of Raouafi & Solanki (2004, 2006) appear to dif- 
fer from earlier empirical modeling efforts, we constructed 
two main sets of electron and flux tube parameters: model R, 
which is designed to replicate many of the conditions assumed 
by Raouafi & Solanki (2004, 2006), and model C, which is es- 
sentially the same as used by Cranmer et al. (1999). Below, 
we also discuss hybrid models with various combinations of 
the conditions assumed in models R and C. 

Model C uses an electron temperature derived by Ko et al. 
(1997) from measurements of ion charge states in the fast so- 

2 Because the degree of radial increase in Tjj_ is relatively uncertain, we 
constructed an additional set of empirical models with no radial increase in 
7jj_ (i.e., where the plane-of-sky values were kept constant over the LOS). 
The resulting probability distributions (§ 4) were virtually identical to those 
computed with the specified radial increase along the LOS. 
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lar wind made by the SWICS instrument on Ulysses (Gloeck- 
ler et al. 1992). We utilize the fitting formula 



T e (r) 



10 6 K 



035 % 



i.i 



+ 1.9 — 



-6.6' 



-i -1 



(5) 



For the electron density, model C uses the expression de- 
rived by Cranmer et al. (1999) from direct inversion of 
UVCS/SOHO white-light polarization brightness (pB) data 
over the poles at solar minimum; i.e., 



nAr) = 3890 ^ 
10 5 cm- 3 V r 



10.5 



+ 8.69 1 ^ 



2.57 



(6) 



The above electron density is a mean value for polar coro- 
nal holes (intermediate between plumes and interplume re- 
gions) between r w 1.5 and 4 Rq. Model C also uses the 
three-parameter empirical function of Kopp & Holzer (1976) 
to specify the superradial expansion of a polar coronal hole. 
The transverse area A(r) oc r 2 f(r) of the entire coronal hole is 
specified by 



/W = l + (/max-l) 



l-exp[(/? -r)/<7i] 
l+exp[(/?!-r)/CTi] 



(7) 



and Cranmer et al. (1999) determined / max = 6.5, R\ = 1.5R®, 
and <j\ = 0.6Rq. Also, the area of the hole is normalized by 
setting the basal colatitude 0o to 28°. The field lines inside 
the coronal hole volume are assumed to self-similarly follow 
colatitudes that remain proportional to the overall boundary of 
the coronal hole at any given radius (see Cranmer et al. 1999). 

Model R uses a constant electron temperature of 10 6 K. 
This value is lower than the peak of the Ko et al. (1997) model 
(T e w 1.5 x 10 6 K at r w 1.6fl Q ), and higher than the value 
from this model at the coronal base (T e w 4 x 10 5 K at r = Rq). 
The constant value of 10 6 K seems to be in closer agreement 
to both theoretical models that take account of strong electron 
heat conduction in the corona (e.g., Lie-Svendsen & Esser 
2005; Cranmer et al. 2007) and with SUMERI SOHO obser- 
vations made above the limb in coronal holes (e.g., Wilhelm 
et al. 1998; Doschek et al. 2001). 3 For the electron density, 
model R uses equation (2) of Doyle et al. (1999), i.e., 



njf) 
10 5 cnr 



= 1000 



+ 0.025 [ ^ 



+ 2.9 



(8) 

To specify the superradial geometry of flux tubes in the polar 
coronal hole, model R uses the analytic magnetic field model 
of Banaszkiewicz et al. (1998). 

Note that Raouafi & Solanki (2004, 2006) used equation (1) 
of Doyle et al. (1999), which is a one-parameter hydrostatic 
fit to various measured electron densities. Above a height of 
r w 6Rq, though, the radial decrease of n e in the hydrostatic 
expression becomes substantially shallower than an inverse- 
square radial decrease. This is not generally expected to oc- 
cur; i.e., in most observations and models, the radial decrease 
in n e goes from a rate much steeper than l/r 2 at low heights 
to l/r 2 at large heights where the geometry is radial and the 
wind speed is constant. The shallow radial density decrease in 
a hydrostatic model is probably unphysical and could lead to 

3 The discrepancies between electron temperatures derived from spec- 
troscopy and from frozen-in ion charge states are not yet fully understood 
(see, e.g., Esser & Edgar 2000, 2001; Chen et al. 2003; Laming & Lepri 
2007). 
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FIG. 2. — Magnetic field lines in the plane of the sky for (a) the B98 (Ba- 
naszkiewicz et al. 1998) model, and (b) the C99 (Cranmer et al. 1999) model 
that used the Kopp & Holzer (1976) flux-tube area function. The dotted hor- 
izontal line denotes the position of the LOS along which various quantities 
are plotted in (c). In panel (c), the superradial angle S is given as a function 
of |;c| (it is the same in the foreground and background halves of the LOS) for 
the two models shown above (see labels). 

an overestimated contribution from large distances along the 
LOS. 

Figure 2 illustrates the differences between the magnetic 
geometries used in models R and C. Figures 2a and 2b show 
field lines that are distributed evenly in polar angle 9 between 
0° and 29° as measured on the solar surface from the north 
pole. The superradial angle 5 characterizes the departure from 
the radial direction, and it is shown in Figure 2c as a function 
of LOS distance x for a polar observing height z = 2.5Rq. 
Formally, 6 is defined as the angle between the radius vector r 
and the magnetic field B (assuming the field points outward), 
i.e., 

'"^(fflSi)- <9) 

For the polar observations described here, the LOS projec- 
tion of any quantity that follows the magnetic field (e.g., the 
outflow velocity) is given by multiplying its magnitude by 
sin(# + <5). The Banaszkiewicz et al. (1998) model exhibits a 
larger degree of departure from radial geometry than does the 
Cranmer et al. (1999) model. However, at the large heights 
for which the UVCS O VI anisotropy results are of interest 
here, the relative differences between the two models — and 
also the differences between either model and a radial geom- 
etry (5 = 0) — are small. 

Figure 3 shows the range of electron densities measured by 
several instruments in polar coronal holes. The strong radial 
decrease in n e (r) has been removed by dividing all measure- 
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FIG. 3. — Comparison of measured electron densities in polar coronal 
holes. Values of n e from Guhathakurta & Holzer (1994) (thick dashed line), 
Fisher & Guhathakurta (1995) (dotted line & vertical bar), Guhathakurta et 
al. (1999) (dot-dashed line & vertical bar), and Cranmer et al. (1999) (thick 
solid line) were divided by equation {8), i.e., equation (2) of Doyle et al. 
(1999). The polar theoretical model of Cranmer et al. (2007) is also shown 
(thin dashed line), as is the approximate hydrostatic fit from eq. (1) of Doyle 
et al. (1999) (triple-dot-dashed line). Gray-scale histogram boxes show the 
range of n e values from the plume statistics study of Cranmer et al. (1999), 
with darker shades denoting the most likely values at each height. 

ments by equation ©. The use of this normalization more 
clearly illustrates the relative differences between the differ- 
ent sets of values, which Raouafi & Solanki (2004, 2006) 
claimed to be important in the derivation of 5+ temperature 
anisotropy. The differences between plumes and interplume 
regions is certainly responsible for some of the wide range of 
variation, but some of it may also be due to absolute calibra- 
tion uncertainties between instruments. Note, though, that the 
curve representing the hydrostatic equation (1) of Doyle et al. 
(1999) appears to be clearly divergent from the other empiri- 
cal curves above r w 8P©, with a slope that is flatter than the 
other measurements even several solar radii below that. 

The gray-scale histogram boxes in Figure 3 illustrate the 
variations due to differing concentrations of polar plumes 
along a single polar LOS over three months (1996 Novem- 
ber 1 to 1997 February 1) using a consistent data set and large 
enough count rates to make Poisson uncertainties negligible 
(see Table 3 of Cranmer et al. 1999). The curves from Fisher 
& Guhathakurta (1995) and Guhathakurta et al. (1999) show 
averages of the various plume and interplume values given in 
those papers, with vertical lines illustrating the relative con- 
trast between the densest plume-filled lines of sight and the 
regions with the fewest numbers of plumes. (These vertical 
lines are shown at r rs 10/?© for clarity, but they are represen- 
tative of the values at the lower heights corresponding to the 
observed white-light data.) Polar electron density values re- 
ported recently by Quemerais et al. (2007) are not shown, but 
they are similar in radial shape to the Fisher & Guhathakurta 
(1995) mean curve (but with values about 10%-20% higher). 
Overall, the variations between data sets that appear to exceed 
the plume-interplume contrast may be due to different instru- 
ment calibrations. 

When modeling the O VI observations summarized in § 2, 
it is probably incorrect to use the lowest "pure interplume" 



electron densities. At z ~ 2.5-3 Rq in polar coronal holes, 
the UVCS observations were typically integrated over 15' to 
30' in the tangential direction, whereas polar plumes at these 
heights have transverse sizes of only about 1' to 2'. Thus, 
the most appropriate electron densities to use are those that 
average over plumes and interplume regions. The lower lim- 
its from Fisher & Guhathakurta (1995) and Guhathakurta et 
al. (1999), as well as the fitting function given by Esser et 
al. (1999), seem to be inappropriate to apply to the empiri- 
cal modeling of these UVCS O VI data. At lower heights, 
where the plume and interplume regions have been resolved 
by UVCS (e.g., Kohl et al. 1997; Giordano et al. 2000), the 
use of the full range of plume and interplume values of n e 
would be warranted. 

3.3. Comparison with Observations 

Once a model data cube (which varies un\, 7]j_, and 7}x/7;|| 
along its axes) has been produced for a given observing height 
z and a given set of n e , T e , and flux tube parameters, the next 
step is to compute the probability of agreement between a 
given observation and each of the simulated observations in 
the cube. We compute this probability P as the product of 
two quantities that are assumed to be independent of one an- 
other: ( 1 ) the probability Pr that the observed line ratio agrees 
with the simulated ratio, and (2) the probability Ps that the 
observed profile shape of the O VI A 1032 line agrees with the 
simulated shape. Because the brighter O VI A1032 line tends 
to dominate the measured "weighted" line width V\/ e , we use 
only the simulated O VI A 1032 line shape in the latter com- 
parison. 

The line ratio probability Pr is relatively straightforward to 
compute. The modeled total intensities of the two compo- 
nents of the doublet are determined by summing up the spe- 
cific intensities over the 12 wavelength bins. Their ratio thus 
gives 7?. mo dei. The relative distance between 7?. mo d e i and the 
observed ratio 7?. b s , in units of the observational standard de- 
viation {SIZohs), is the quantity that determines the probability 
of agreement. Assuming the uncertainties are normally dis- 
tributed, the probability is 



1-erf 



_ 'rZ: 



model | 



STZobs \/2 



(10) 



(see, e.g., Bevington & Robinson 2003). A larger argument in 
the error function ("erf" above) denotes a larger discrepancy 
between the modeled and observed ratios, and thus a lower 
probability of agreement. 

The line shape probability Ps is not as easy to compute as 
Pr. An initial attempt was made to fit the simulated profiles 
with Gaussian functions, and then to compare the resulting 
V\i e widths with the observed values using a similar expres- 
sion as equation ( [10) . However, there were many instances 
where the modeled lines were far from Gaussian in shape, 
but the best-fitting Gaussian (which was a poor fit in an ab- 
solute sense) happened to agree with the observed Vj i e . This 
resulted in spuriously high probabilities for wide regions of 
parameter space that should have been excluded. Thus, we 
found that the tabulated specific intensities (i.e., the full line 
shapes) need to be compared on a wavelength-by- wavelength 
basis. This raises the issue of what to use for the "observed" 
line shape. As described in § 2, the UVCS/SOHO data points 
contain a wide range of instrumental effects that were taken 
into account in the line fitting process. In order to compare 
similar quantities, either these effects must be inserted into the 
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model profiles, or we must reconstruct "observed profile" in- 
formation from the extracted V\/ e measurements and the SV\/ e 
uncertainties. We chose the latter option. 

To determine the probability of agreement between the set 
of modeled specific intensities (/A.modei) and the reconstructed 
observed intensities (/A,obs)> we computed a \ 2 quantity, 

1 -1 n2 

^A.obs 1 X, model 



x 



= E 



Sh 



,obs 



(ID 



where /A.modei came from the data cube, and /A.obs was con- 
strained to be a Gaussian function with the observed V\/ e 
width and a total intensity equal to that of the modeled profile. 
(The observed total intensity was not used because the com- 
parison being done here is only between the relative shapes.) 
The <5/A,obs uncertainty was computed as a function of wave- 
length by comparing the idealized /A.obs profile with two oth- 
ers computed with line widths Vy e -5Vi/ e and Vi/ e + 5Vi/ e 
(with all three profiles normalized to the same modeled to- 
tal intensity). These three profiles exhibited a range of spe- 
cific intensities at each wavelength, and the standard deviation 
quantity Sl\ obs was defined as half of that full range. Then, 
the x 2 quantity above constrains the probability that the ob- 
served and modeled profiles are in agreement (i.e., the proba- 
bility that the observed and modeled specific intensity values 
are drawn from the same distribution). Assuming normally 
distributed uncertainties, this probability is given by 

1 



= Q(X» = 



dt 



(12) 



I>/2), xV2 

(Press et al. 1992), where v = N\ - 1 is the effective num- 
ber of degrees of freedom (for N\ = 12 wavelength points) 
and T(x) is the complete Gamma function. When x 2 <^ v the 
above probability approaches unity (i.e., the modeled profile 
is a good match to the observed profile), and when x 2 3> v the 
above probability is negligibly small. 

We thus obtained the total probability P = PrPs as a func- 
tion of the three main 5+ variables of each data cube, for each 
observation at the height z consistent with that data cube. The 
question of what is considered to be a large or small proba- 
bility is open to some interpretation. Below, we often use a 
standard "one sigma' ' probability P Xa = 1 -erf(l /\/2) = 0.317 
as a fiducial value above which the solutions are considered 
to be good matches with the data. 

4. EMPIRICAL MODEL RESULTS 

In this section we present results for the most probable val- 
ues of the 5+ ion properties between 1.5 and 3.5 Rq. In 
§ 4.1, we show how the essential information inside the three- 
dimensional probability cubes can be extracted and analyzed 
in a manageable way. In § 4.2, the optimal values for O s+ out- 
flow speed and perpendicular temperature are presented for 
models C and R. The resulting values of Uju and 7/j_ are con- 
sistent with earlier determinations of preferential ion heating 
and acceleration with respect to protons. In § 4.3, we discuss 
the determination of the anisotropy ratio T t ± /T^ for models C 
and R, which is less certain than the other two quantities. We 
then focus in detail on a single representative height in § 4.4 
in order to determine how models C and R can give rise to 
qualitatively different conclusions about the ion temperature 
anisotropy. Finally, in § 4.5 we extract information from both 
models C and R about the 5+ ion concentration in the ex- 
tended corona — i.e., we compute the ratio n s+/nu from the 
comparison of observed and modeled O VI A 1032 total inten- 
sities. 
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FIG. 4. — Contour plots of the maximum probabilities of agreement be- 
tween the model R data cube and the UVCS observation from 1997 January 
5 at r = 3.07Sq. The three panels show probabilities as a function of each 
unique pair of the three data-cube axis quantities. Contour levels are plotted 
at 90% of the probability values 1, 0.3, 0.1, 0.03, and 0.01 (see gray-scale 
code in panel a). 

4.1. Deriving Ion Properties from the Data Cubes 

We constructed two sets of radially dependent data cubes: 
one for the model C assumptions for n e , T e , and flux-tube ge- 
ometry, and one for model R. Each set consisted of 13 data 
cubes with observation heights z= 1.5, 1.6, 1.7, 1.8, 1.98, 
2.11, 2.3, 2.42, 2.563, 2.7, 3.0, 3.09, and 3.565 R e . These 
values were chosen to align with the observed data points 
shown in Figure 1 . Any discrepancies between the observed 
and modeled heights never exceeded ±0.065 Rq, and for the 
whole data set the average absolute value of the discrepancy 
was only 0.012 Rq. We then constructed 49 individual "prob- 
ability cubes" for each of the data points for which both V\ / e 
and 1Z exist. 

Even for just a single comparison between an observation 
and a data cube at one height, it is a challenge to display 
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the full three-dimensional nature of the probability "cloud" 
P(W;|i,T;j_,7}j_/7'i||). We limit ourselves to showing lower- 
dimensional projections that keep only the highest probabil- 
ity values taken over the axes that are not being shown. As 
an example, in Figure 4 we display two-dimensional contours 
of P as a function of all three unique pairings of the three 
axis-quantities of the data cube. The specific comparison is 
between the measurement shown in Table 1 from 1997 Jan- 
uary 5 (3.O7/? , Vy e = 690 km s" 1 ) and the model R data 
cube constructed at z = 3.09Pq. In all three contour plots, the 
probability shown at each location is a maximum taken over 
the third quantity that is orthogonal to the projection plane. 
Thus, for regions with a low probability in these diagrams, 
we can be assured that there are no values of the unplotted 
coordinate that can give synthetic line profiles in agreement 
with the observations. 

Figure 4a shows an approximate anticorrelation between 
the ion outflow speed and the perpendicular kinetic temper- 
ature in the subset of generally "successful" models. This 
arises mainly because the lines can be broadened both by mi- 
croscopic LOS motions (roughly proportional to 7}j_) and by 
the projection of the superradially flowing bulk outflow speed 
along the LOS (which goes as Uj\\), When one of these quanti- 
ties goes up, the other must go down in order to match a given 
observed line profile. Figure 4b shows that the region of pa- 
rameter space with the larger contribution by Ti± (upper left) 
also requires a large anisotropy ratio, but the region with the 
larger contribution by bulk LOS motions (lower right) may 
be able to match the observations with an isotropic velocity 
distribution (i.e., 7}j_/7/|| w 1). 

The large amount of information in contour plots like Fig- 
ure 4 can be collapsed down to a smaller list of parame- 
ters. We created three one-dimensional probability curves as 
a function of each of the three main axis quantities, with the 
maximum values extracted from the full plane subtended by 
the remaining two neglected quantities. Thus, we define the 
reduced probability functions P„(w,|| ), P< and P a (7}j_ /7/m) 
(the subscript "a" denotes the anisotropy ratio). These func- 
tions are generally peaked at some high value close to 1 and 
exhibit lower values far from the optimal solutions. Figure 5 
shows these reduced probabilities for the same example case 
shown in Figure 4. The reduced probabilities for Uf m and 7}j_ 
are peaked relatively sharply around their most probable val- 
ues. Note that we plot the perpendicular kinetic temperature 
in units of a most-probable speed Wi± = QkgJix/md 1 ' 2 m or- 
der to facilitate comparison with earlier papers. The reduced 
probability for the anisotropy ratio, shown in Figure 5b, is 
less centrally peaked and thus the best solution for this value 
is less certain. The peak value corresponds to a most-probable 
anisotropy ratio of 7jj_/7fii w 6, but note that the probability 
of isotropy remains reasonably high at ^40%. 

It is interesting to contrast the exhaustive data-cube-search 
technique used here with the more straightforward approaches 
taken in earlier papers. For example, Raouafi & Solanki 
(2004, 2006) simulated the properties of the O VI AA1032, 
1037 lines after first fixing the radial variation of the outflow 
speed and ion temperature. Figure 5b shows an illustrative 
"cut" through the data cube at fixed values of uqi = 600 km s" 1 
and Wi± =215 km s -1 at r = 3.097? Q (similar to the values 
used by Raouafi & Solanki 2006). In this case, the most prob- 
able value of the anisotropy ratio is surprisingly close to 1, as 
was also assumed by Raouafi & Solanki. The apparent consis- 
tency with the observations (i.e., a value of P a of about 30%) 
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FIG. 5. — Reduced probabilities for one specific comparison between a 
UVCS observation at r = 3.07 Rq and the model R data cube (see also Figs. 4, 
8, 10, and 11). (a) P„ versus outflow speed Ufa (solid line) and P t versus per- 
pendicular most-probable speed w,x (dashed line), (b) P a versus anisotropy 
ratio Tj±/T}\\ for a search of the entire data cube (solid line) and for a "slice" 
through the data cube with fixed values of and Wj± given above (dashed 
line). Also shown is the threshold level P\ a (dotted lines) defined in the text. 

may be misleading if the rest of the data cube is not searched. 
Thus, we can assert that any results concerning the anisotropy 
ratio that were obtained by not searching the full range of pos- 
sibilities for the ion parameters are potentially inaccurate. 

The ultimate goal of the empirical modeling process is to 
characterize the peak values and widths of the reduced prob- 
ability curves, in order to obtain the optimal measured values 
(with uncertainty limits) for the relevant 5+ plasma proper- 
ties. The most satisfactory outcome, of course, would be very 
narrow peaks that occur far from the edges of the parameter 
space, but this is not always the case. After some experimen- 
tation, we chose to use weighted means to obtain the peak 
values, i.e., 

J dxxPJx) 



(x) = 



(13) 



J dxP x (x) 

where x denotes any of the three axis quantities um, T/±, or 
log(Iix/7i|i). We used the logarithm of the anisotropy ratio 
in equation ( fT3l ) because tests showed that if the ratio itself 
(which spans three orders of magnitude) was used, the above 
mean would be weighted strongly toward the largest values 
even when the maximum of the probability distribution is at 
much lower values. 

We experimented with using the variance, or second mo- 
ment, of the reduced probability distributions to characterize 
the widths of "error bars" for each measurement. However, 
because the probability curves are generally not symmetric 
around the peak values, the second moment often did not 
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FIG. 6. — Derived outflow speeds (a), perpendicular most-probable speeds 
(b), and kinetic anisotropy ratios (c) for model R (red points) and model C 
(blue points). Symbols show the weighted means of the reduced probability 
distributions, with styles the same as in Figure 1. Vertical bars show the 
full range of parameter space with reduced probabilities greater than P\ CT (for 
model R). Also shown are empirical models Bl and B2 from Cranmer et al. 

(1999) (dotted lines) and Alfven wave quantities (Sv^)^ 2 in panel (b) and 
A c ff in panel (c), derived from the model of Cranmer et al. (2007) (dot-dashed 
lines). 

accurately give a range of values with reasonable probabili- 
ties. Instead, we performed a straightforward search for the 
range of probabilities that are higher than the threshold value 
P\ a ~ 0.317 discussed above. The lower and upper limits of 
that range were taken to be the ends of the uncertainty bounds 
for each measurement. 

4.2. Preferential Ion Heating and Acceleration 

Figure 6 shows the weighted mean and error-bar quantities 
for the O s+ plasma properties, defined as in the previous sub- 



section, as a function of heliocentric height. The results from 
model C and model R are plotted in two different colors, with 
only the error bars of model R shown for clarity. Here we fo- 
cus on the ion outflow speed (Fig. 6a) and the perpendicular 
ion temperature (Fig. 6b). On average, the derived values of 
(k;||) and (w;_l) were consistent between the two sets of mod- 
els. To quantify the impact of varying the electron density, 
electron temperature, and flux-tube geometry, we computed 
ratios of the model C values to the model R values for each 
data point. For the 49 data points taken together, the mean 
value of the ratio of the outflow speeds was 1 .002, with a stan- 
dard deviation of 19%, and the mean value of the ratio of per- 
pendicular most-probable speeds was 0.991, with a standard 
deviation of 14%. This shows that the determination of these 
parameters is relatively insensitive to the choices of electron 
density, electron temperature, and flux-tube geometry. 

The radial dependence of the derived (Uj») values in Fig- 
ure 6a is similar to that of the O s+ empirical models B 1 and 
B2 given by Cranmer et al. (1999). Note the emergence of a 
natural trend of radial acceleration in (w;ii), with the possible 
exception of the data points at r > 3.5Rq. This is especially 
serendipitous given that each data point was analyzed inde- 
pendently of all others. 

The derived 5+ outflow speeds support earlier claims of 
preferential ion acceleration in coronal holes. At r = 2.5Rq, 
the range of ion outflow speeds that gives rise to high prob- 
abilities of agreement with the data points is approximately 
280-500 km s" 1 At this height, these values are substantially 
larger than bulk (proton-electron) solar wind outflow speeds 
derived via mass flux conservation. Figure 41 of Kohl et al. 
(2006) showed a selection of 12 bulk outflow speed models 
derived using all possible combinations of four n e models and 
three coronal-hole geometries. At r = 2.5 Rq, these 12 models 
gave a range of bulk outflow speeds of 1 15-300 km s" 1 . De- 
spite the small degree of overlap between the two ranges, the 
mean value of the O s+ range (390 km s -1 ) exceeds the mean 
value of the mass flux conservation range (208 km s" 1 ) by 
almost a factor of two. Also, proton outflow speeds derived 
from H I Lya Doppler dimming (from a selection of papers 
all dealing with polar coronal holes at the 1996-1997 solar 
minimum) were shown in Figure 41 of Kohl et al. (2006). At 
2.5 Rq, the range of these values is 160-260 km s" 1 — with a 
mean of 210 km s -1 — which is still significantly lower than 
the range of O s+ ion outflow speeds discussed above. 

In Figure 6b, the trend of radial increase in (w;j_) is also 
roughly similar to that found by Cranmer et al. (1999), espe- 
cially below about 2.3 Rq. At larger heights, though, there 
appears to be less evidence for a systematic increase than ex- 
isted in the model Bl and B2 curves. This could have re- 
sulted either from the inclusion of the new data points or from 
the more exhaustive treatment of uncertainties in the new pa- 
rameter determination method described above. However, if 
one takes all of the derived (w,jj values for r > 2Rq and fits 
to a straight line, the best-fitting slope is still increasing with 
height at a rate of 50 km s" 1 per Rq. This is about a third of 
the ~150 km s" 1 per Rq slope in the Bl and B2 models. 

The most-probable speeds shown in Figure 6b, al- 

though slightly smaller than those given by the Cranmer et 
al. (1999) model Bl and B2 curves at some heights, still show 
definite evidence for preferential ion heating. The mean value 
of the (vfjx) values at heights r > 2.5Rq in Figure 6b is 363 
km s" 1 , with a standard deviation of 73 km s" 1 . Between 2.5 
and 3 Rq, the perpendicular proton most-probable speeds de- 
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rived from H I Lya were about 2 10-240 km s -1 , with a mean 
value of about 225 km s -1 (see models Al and A2 of Cranmer 
et al. 1999). The fact that the O s+ mean value exceeds the 
proton mean value by almost two standard deviations implies 
that the 5+ kinetic temperature at this height is very likely 
to be more than "mass proportional" (i.e., implying an oxy- 
gen kinetic temperature of 130 MK, or more than 40 times 
the proton kinetic temperature of ~3 MK. 

It is important to note that the derived kinetic temperatures 
are likely to be a combination of thermal and nonthermal mo- 
tions. The ion-to-proton kinetic temperature ratio of ^40, 
derived above, is likely to be a lower limit to the true ratio 
of thermal, or microscopic temperatures. If unresolved wave 
motions are deconvolved from the empirical values of (w,-j_), 
the proton most-probable speed will be reduced by a relatively 
larger amount than the 5+ speed. As an example, the the- 
oretical polar coronal hole model of Cranmer et al. (2007) 
has a LOS-projected Alfven wave amplitude at r = 2.5 Rq of 
116 km s" 1 (this is also plotted in Fig. 6b). Converting these 
motions into temperature-like units and subtracting from both 
values given above, one obtains an 5+ perpendicular temper- 
ature of 115 MK and a proton perpendicular temperature of 
2.2 MK. The ratio of ion to proton temperatures has thus in- 
creased from about 40 to 50. In any case, it is clear that the 
dominant contributor to the ion kinetic temperature is the true 
"thermal" temperature, with only a relatively minor impact 
from broadening due to macroscopic motions. 

4.3. The Ion Anisotropy Ratio 

Figure 6c illustrates the largest discrepancy between the 
empirical models of Cranmer et al. (1999) and the present 
models (both C and R). Above a height of ^2.5 Rq, models 
Bl and B2 demanded a strong anisotropy ratio Ti±/T^ > 10, 
but the optimal ratios derived in this paper seem to cluster 
between 2 and 10 with no discernible radial dependence. It 
is important to note, though, that there is considerable over- 
lap of the uncertainties between the old and new ranges of 
(7}j_/7]||). Several of the error bars shown in Figure 6c extend 
up into the range of ratios from models Bl and B2. Also, the 
dotted curves that illustrate models B 1 and B2 correspond to 
the "optimal" values of the anisotropy ratio from Kohl et al. 
(1998) and Cranmer et al. (1999); the uncertainties in those 
models are not shown. 

Because the derived values of the kinetic temperature ratio 
(7;j_/7]|i) in Figure 6c exceed unity by only a relatively small 
amount, it is worthwhile to examine whether the numerator 
(Ti±) may have been enhanced by unresolved wave motions 
perpendicular to the magnetic field. In other words, for a real- 
istic model of perpendicular wave amplitudes in polar coronal 
holes, we investigate whether a truly isotropic microscopic 
velocity distribution could have given an effective anisotropy 
ratio that exceeds 1 . We compute such an effective anisotropy 
ratio as 

Aeff = 1 771 2 \ / 2 x ' ( 14 ) 

l-((Svi) x /wf ± ) 

where (Sv 2 ± ) x is the square of the frequency-integrated Alfven 
wave amplitude divided by two to sample only the motions 
along one of the two perpendicular directions (i.e., only along 
the LOS or x axis). As above, we used the Alfven wave prop- 
erties from the turbulence-driven polar coronal hole model of 
Cranmer et al. (2007). The model wave amplitude is plotted 
in Figure 6b and the quantity A e g is plotted in Figure 6c. The 
condition A e g w 1 corresponds to the situation where the am- 
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FIG. 7. — Probability of ion isotropy P a (l) plotted versus heliocentric dis- 
tance for model R, using the same data symbols as Figures 1 and 6. The filled 
circle shows the probability of isotropy for the specific observation (3.07 Rq, 



= 690 km s ') that is considered in detail in Figures 4, 5, 8, 10, and 11. 



plitudes are too small to contribute to the anisotropy ratio (as 
defined in the empirical models). Below rsl .7 Rq, the curve 
in Figure 6c shows that A e s does indeed exceed 1 by about the 
amount computed from the UVCS data. At these low heights, 
the derived value of (w,-j_) is of the same order of magnitude 
as the wave amplitude, so the latter can "contaminate" the 
determination of the true perpendicular most-probable speed. 
At heights larger than about 2 Rq, though, the wave ampli- 
tudes are small in comparison to the derived (w,-j_) values, and 
thus A e ff w 1. We thus conclude that above 2 Rq, any derived 
anisotropy ratio (Tj±/T^) is likely to be truly representative 
of the microscopic velocity distribution and not affected by 
wave motions. 

Despite the comparatively low values of the anisotropy ra- 
tio shown in Figure 6c ((7j-j_/7/||) w 2-10), we should em- 
phasize that these values are often significantly different from 
unity. It is important to note that all of the data points have 
their largest reduced probability — measured either using the 
weighted mean defined above or by simply locating the max- 
imum value — for anisotropy ratios larger than unity. 

The preponderance of evidence for anisotropy is also il- 
lustrated in Figure 7, which shows the probability that each 
measurement could be explained by an isotropic 5+ veloc- 
ity distribution. In other words, Figure 7 gives the value of 
P a (l) for each probability cube. Taken together, a signifi- 
cant majority of the values (about 78% of the total number) 
fall below the fiducial one-sigma value of P\ a , indicating that 
isotropy should not be considered a "baseline" assumption. 
Below about r = 2Rq, a few of the measurements correspond 
to large probabilities that an isotropic distribution can explain 
the observations. Note from Figure 6c, though, that the most- 
probable anisotropy ratios for these measurements tend to be 
greater than 1, but some of the error bars extend down past 
T t ± /7}|| = 1 . However, between 2. 1 and 2.7 Rq the probability 
of isotropy is very small for all of the observed data points. 
Above 3 Rq, some of the values of P a (l) become large again, 
but we believe this may be due to the relatively high observa- 
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FIG. 8. — Reduced probability P a versus ion anisotropy ratio 7}x/7]|| for 
one specific comparison between a UVCS observation at r = 3.07 Rq and the 
model R data cube (see also Figs. 4, 10, and 11). The probability computed 
with the actual observational uncertainties (thick solid line) is compared with 
trial curves computed with a range of constant factors e multiplying 8V l i e 
and S7Z (dashed lines). Also shown is the threshold level P\ a (dotted line) 
and the probability of isotropy for the standard e = 1 case that is also shown 
in Figure 7 (filled circle). 

tional uncertainties on the O VI intensities and line widths at 
these large heights (see below). 

To better understand the impact of observational uncer- 
tainties on the probability of isotropy, Figure 8 shows the 
full P a (7ij_/7]-||) curves for one specific measurement at r = 
3.07/? Q (i.e., the same measurement used in Fig. 5). The 
multiple curves were constructed by multiplying the known 
observational uncertainties 6V\/ e and 51Z by arbitrary factors 
e. Generally, larger uncertainties lead to lower \ 2 values when 
comparing the observed and modeled line shapes, and thus to 
larger probabilities of agreement between the observed and 
modeled profiles. Interestingly, though, the anisotropy ratio 
Ti±/TjH at which the maximum probability occurs remains 
roughly constant when e is varied between 0.5 and 2. Thus, 
if future observations above 3 Rq were to obtain the same 
general range of values for V\/ e and 1Z but with lower uncer- 
tainties, it could provide stronger evidence for ion anisotropy 
up at these heights. 

Although Figure 6c does not seem to indicate a substantial 
difference between models R and C, it is useful to compare 
these models in some additional detail. For all data points, 
the mean ratio of model C to model R anisotropy ratios was 
1.199, but the large standard deviation (76%) shows that the 
models are often quite different from one another. Taking only 
the heights above 2.2 R & , the mean ratio of model C to model 
R anisotropy ratios increases to 1.514, indicating that on aver- 
age model C generates larger anisotropics than model R over 
the height range where anisotropies appear to be required. 

Figure 9 shows the ratio of model C to model R anisotropy 
ratios as a function of height. Below r w 2.2R Q the two 
models produce roughly the same result for the anisotropy ra- 
tio. Above that height, the solutions split into two groups: 
one where model C produces a substantially larger ratio (2- 
3 times that of model R), and one where model C produces 
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FIG. 9. — Ratios of model C to model R values for the weighted mean 
anisotropy ratios (7}x/Zj|| ) shown as a function of height and using the same 
data symbols as Figures 1, 6, and 7. 

a comparable or slightly smaller ratio than model R. Note 
that the height range of 3.0-3.1 Rq — over which model R 
predicted a rise in the probability of isotropy (see Fig. 7) — 
strongly favors larger anisotropies for model C. 

4.4. Varying the Electron Density, Electron Temperature, and 

Geometry 

One of the main motivations for this paper was to explore 
why the results of Raouafi & Solanki (2004, 2006) were so 
different from earlier results (e.g., Cranmer et al. 1999) re- 
garding the 5+ anisotropy ratio. In this subsection, we study 
the differences between model R and model C in more de- 
tail by focusing on the shapes of the reduced probability dis- 
tributions for one representative data point. As in Figures 
4, 5, and 8, we used the probabilities generated by com- 
paring the WCS/SOHO measurement from 1997 January 5 
(r = 3.07/?q, V\i e = 690 km s" 1 ) with data cubes constructed 
with various assumptions. This data point is denoted by a 
filled circle in Figure 9, and it is clear that this point is rep- 
resentative of the majority of the data points (5 out of 7) at 
r^3R Q . 

Figure 10 shows a range of reduced probability curves 
Pa(Ti±/Tjii) that were computed from data cubes constructed 
with various combinations of the model R and model C pa- 
rameters. The three-letter names for the models denote the 
individual choices for n e , flux-tube geometry, and T e (in that 
order). The "pure" model R and model C cases are thus called 
RRR and CCC. 

Before examining the impact of the individual parameters 
on the reduced probability curves, we note that the model 
CCC curve in Figure \0b mirrors almost exactly the results 
of Kohl et al. (1998) and Cranmer et al. (1999) at r w 3fl : 
the most likely 5+ anisotropy ratio ranges between 10 and 
100, and an isotropic distribution is highly improbable. Model 
CCC exhibits a most probable ion outflow speed (k,-m) = 508 
km s" 1 , which is only marginally smaller than the model RRR 
value of 521 km s" 1 . Model CCC has an optimal solution for 
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FIG. 10. — Reduced probability P a versus ion anisotropy ratio 7jj_/7}|i for 
the same data comparison as in Figures 4 and 8, but for various combinations 
of the model R and model C parameters (see above for line styles). The 
order of the three-letter designations is {n e ,f(r), T e }. Panel (a) thus shows 
all models computed with the Doyle et al. (1999) n e and panel (b) shows 
all models computed with the Cranmer et al. (1999) n e . Also shown is the 
threshold level P\,y (horizontal dotted lines) and the probability of isotropy 
for the main RRR and CCC cases (filled circles). 

(wi±), though, of 541 km s" 1 , which is 23% larger than the 
corresponding value of 440 km s" 1 for model RRR (i.e., a 51% 
higher value of (T i± _) for model CCC). Model CCC tended to 
produce more line broadening via "thermal" motions near the 
plane of the sky, and model RRR tended to produce more line 
broadening via bulk outflow projected along the LOS. 

The other curves shown in Figure 10 explore which of the 
three varied parameters were most responsible for the differ- 
ences between models RRR and CCC. We see immediately 
that the choice of electron temperature T e , which in our mod- 
els impacts only the collision rate q\2, is relatively unimpor- 
tant. The 8 curves can thus be separated into 4 pairs, each 
of which has the same choice for n e and flux-tube geometry 
(i.e., RRX, RCX, CRX, and CCX, where 'X' denotes either 
option for 7^). The overall insensitivity to electron tempera- 
ture is evident from the fact that the two curves in each pair 
are virtually indistinguishable from one another. 

Figure 10 shows that the unique features of the RRX mod- 
els (i.e., a higher probability of isotropy and a strong peak at 
Tix/Tj\\ < 10) are only present when both the electron density 
and flux-tube geometry are treated using model R. The mod- 
els with only one of these two parameters treated using model 
R (i.e., RCX and CRX) appear more similar to the CCX mod- 
els. At large values of the anisotropy ratio, both the RCX and 
CRX models are virtually identical to the CCX models. At 



low values of the anisotropy ratio, the CRX model is roughly 
intermediate between the CCX and RRX models. Generally, 
though, the combination of the model R assumptions for elec- 
tron density (e.g., Doyle et al. 1999) and flux-tube geometry 
(e.g., Banaszkiewicz et al. 1998) are needed to produce broad 
enough profiles via outflow speed projection along the LOS to 
explain the observations without the need for extreme temper- 
ature anisotropics. Specifically, this enhanced LOS projection 
effect arises for two coupled reasons. 

1. As seen in Figure 2c, the Banaszkiewicz et al. (1998) 
flux tubes are tilted to a greater degree away from the 
radial direction than the Cranmer et al. (1999) flux 
tubes. Because of these larger values of 5, a larger frac- 
tion of the outflow speed u t H is projected into the LOS 
direction (when |x| > 0) for model R. 

2. Figure 3 shows that the Doyle et al. (1999) electron 
density does not drop as rapidly with increasing height 
(between about 3 and 10 Rq) as nearly all of the other 
plotted n e functions. Thus, for observation heights at 
about 3 Rq, the Doyle et al. (model R) electron den- 
sity provides a relative enhancement for points along 
the foreground and background (|jic| > 0) in comparison 
to the plane of the sky (x = 0). 

Note also from Figure 3 that the electron density used for 
model C is about 10% to 30% larger than that used for model 
R at the heights of interest (r w Rq). A higher value of n e 
is expected to result in emission lines that are dominated more 
by the collisional component of the emissivity, which scales 
as n 2 e (eq. |2]), with a correspondingly weaker contribution 
from the radiative component, which scales linearly with n e 
(eq. Q). Because of the different density dependences, the 
collisional component is not extended as far along the LOS 
as the radiative component. Thus, models with higher den- 
sities would be expected to behave more like model C (with 
emission dominated by the plane of the sky), and models with 
lower densities would be expected to behave more like model 
R (with emission extended over a larger swath of the LOS). 

To explore the effects of varying the electron density, we re- 
peated the model R data cube analysis (for the fiducial height 
shown in Figs. 8 and 10) with n e (r) multiplied by constant fac- 
tors. Figure 1 1 shows the resulting reduced probability curves 
as a function of the 5+ temperature anisotropy. A model with 
half of the Doyle et al. (1999) electron density has a lower 
preferred value of T^/Tq and a much higher probability of 
isotropy than the standard model R. A model with double the 
Doyle et al. (1999) electron density resembles model C in that 
there is a high preferred range of Ti±/T t ii and a low probabil- 
ity of isotropy. Despite the large change in appearance of the 
P a curves as shown in Figure 1 1 , the preferred values of the 
outflow speed and perpendicular kinetic temperature do not 
vary by very much as n e is varied up and down by a factor of 
two: changes by only about ±8% (increasing as n e de- 
creases), and (wi±) changes by only about ±10% (increasing 
as n e increases). These determinations appear to be relatively 
insensitive to the choices for n e and flux-tube geometry. 

The ratio of collisional emissivity to the total line emission 
changes dramatically for the models shown in Figure 11. For 
model RRR, the optimal model in the data cube exhibited a 
collisional fraction of 93.7% for the O VI A 1032 line and a 
fraction of 44.6% for the O VI A1037 line (the latter being 
"Doppler pumped"). The model with half of the model R den- 
sity had lower collisional fractions for the AA1032, 1037 lines 
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FIG. 1 1 . — Reduced probability P a versus ion anisotropy ratio 7/x/7]|| for 
the same data comparison as in Figures 4, 5, 8, and 10, but for a range of 
constant multipliers to the Doyle et al. (1999) electron density. The basic 
"model R" n e (solid line) is compared to a model with half {dashed line) and 
double (dot-dashed line) this electron density function. Also shown is the 
threshold level Pi CT (dotted line) and the probabilities of isotropy for the three 
curves (filled and open circles). 



of 88.2% and 28.7%, respectively. The model with double the 
model R density had higher collisional fractions of 96.8% and 
61.7%. 

It is important to note, however, that the differences in colli- 
sionality for the models shown in Figure 10 are not as drastic 
as those shown in Figure 1 1 . Model CCC exhibited collisional 
fractions for the AA1032, 1037 lines of 90.7% and 47.0%. 
These values are only a few percentage points different from 
the model RRR fractions. The other intermediate models have 
values that cluster between those of models RRR and CCC. 
The larger value of n e in the plane of the sky for model C is 
compensated — to some degree — by the slower decrease in n e 
along the LOS for model R. Thus, despite the superficial re- 
semblance between the model CCC curve in Figure 10 and the 
"double n e " curve in Figure 11, one cannot invoke a varying 
amount of collisionality to explain the differences between 
models R and C. The LOS extension effects discussed above 
are more subtle than simply varying n e by a constant amount. 

Another way we explored the dependence of the reduced 
probabilities on electron density was to produce a set of 
three other models with alternate functional forms for n e (r), 
but the same flux-tube geometry and T e as used in model 
R. These models utilized the mean electron density curves 
from Guhathakurta & Holzer (1994), Fisher & Guhathakurta 
(1995), and Guhathakurta et al. (1999) (see also Fig. 3), and 
the O VI data cubes were created only at the fiducial height 
of 3.09 Rq. 4 The reduced probabilities P a for these models 
all fell within the general range of variation illustrated in Fig- 
ure 10 and are not plotted. However, the construction of these 

4 We also created a data cube for the hydrostatic equation (1) of Doyle et 
al. (1999), but this model exhibited an unusually strong extension along the 
LOS. There was a substantial contribution to the O VI emissivity even at the 
LOS integration limits of x = ±15Rq, which actually led to an extremely 
low probability of isotropy. However, we discarded this model because the 
shallow n e at large heights is clearly unphysical. 
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FIG. 12. — Reduced probabilities of isotropy (diamonds) and weighted 
mean anisotropy ratios (triangles) for models having a range of n e values 
and identical flux-tube and T e properties (see text for details). All data-cube 
comparisons were computed for the same fiducial data point illustrated in 
Figures 4, 5, 8, 10, and 11. Curves denote the best fitting quadratic relations 
as a function of the optimized density quantity n\ 83 / n 6' 

models increased the number of data cubes with "model R- 
like" flux-tube and T e parameters to seven: i.e., these three 
new ones, the three models shown in Figure 11, and model 
CRR (with a model C electron density). We performed a 
regression analysis on the seven values of the probability of 
isotropy P a (l) and the weighted mean anisotropy (T^/Tq) to 
find the optimal functional dependence on two "independent" 
variables that characterize the electron density: 



"3 



n e (3R Q ) 



10 6 cm- 3 ' " 6 " 10 6 cm- 3 

where the arbitrary normalizations serve only to keep the 
combined quantities (discussed below) of order unity. The 
quantity H3 characterizes the electron density in the plane 
of the sky of the observation, and the ratio «3 /«6 character- 
izes the large-scale density gradient and thus the relative en- 
hancement of foreground and background regions along the 
LOS. From the discussion above, we expect that larger val- 
ues of both «3 and «3 /«6 should result in lower probabilities 
of isotropy and higher most-probable values of the anisotropy 
ratio. Indeed, the regression analysis found that the modeled 
values of these quantities exhibited the lowest combined \ 2 
spread for a single independent variable that scales as n\ &3 jn^ 
(close to the product of «3 and nj/n^). Figure 12 shows these 
values as well as the best-fitting quadratic functions to P a (l) 
and (7}^/7}||). The combined dependence on both n e and its 
radial gradient appears to be a key factor in determining the 
relative probabilities of isotropy and strong anisotropy. 

Finally, we must evaluate which sets of choices for the elec- 
tron density and the flux-tube geometry are most consistent 
with observed polar coronal holes. Figure 3 does seem to in- 
dicate that most measured n e curves (as well as one exam- 
ple theoretical result for n e ) behave more like the "model C" 
(Cranmer et al. 1999) case than the "model R" (Doyle et al. 
1999) case. Between the heights of about 3 and 10 R Q , the 
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majority of the curves in Figure 3 exhibit a steeper radial de- 
crease than the Doyle et al. (1999) model. Thus, the CCX or 
CRX models shown in Figure \0b appear to be more consis- 
tent with observations than the RRX or RCX models in Fig- 
ure 10a. This then implies that substantial 5+ anisotropy 
(Ti±/Tj\i > 10) is also preferred at large heights. The optimal 
choice of the flux-tube geometry is less certain. Ideally, obser- 
vations of the nonradial shapes of polar plumes should be able 
to constrain the magnetic field geometry (see, e.g., Wang et al. 
2007; Pasachoff et al. 2007), but it is unclear whether exist- 
ing plume observations would be able to distinguish the subtle 
differences between, e.g., Figures 2a and 2b. In any case, the 
geometry does not seem to be as major an issue as the electron 
density, since the variance between the four curves in Figure 
10b is not large. 

4.5. Oxygen Ion Number Density 

By comparing the observed and modeled total intensities 
of the O VI A 1032 line, it is possible to derive firm limits on 
the combined elemental abundance and ionization fraction of 
5+ . The ion concentration is useful both as a tracer of fast 
and slow solar wind streams (e.g., Zurbuchen et al. 2002) and 
as a possible diagnostic of the amount of preferential heat- 
ing deposited in coronal holes (Lie-Svendsen & Esser 2005). 
A first attempt at determining the O s+ number density from 
UVCS data was made by Cranmer et al. (1999), but the "data 
cube search" technique developed in this paper allows a much 
more definitive set of measurements to be made. 

The numerical code that computes the O VI line emission 
used an arbitrary constant value for the ratio /o = n s+/n e 
of 4.959 x 10~ 6 , which was derived from the oxygen abun- 
dance of Anders & Grevesse (1989) and the measured O s+ 
ionization fraction of Wimmer-Schweingruber et al. (1998). 
This is merely a fiducial value that does not affect the fi- 
nal determination of this ratio for a given UVCS observa- 
tion. When comparing the results from an empirical model 
data cube with a specific observation, the probability values 
P(u i u,Ti±,Ti±/T j \\) are used to construct a weighted mean of 
the modeled O VI A 1032 total intensity using equation d!3T >. as 
well as lower and upper limits using the full range of models 
with probabilities that exceed P\ a . These values are converted 
into "observed" ion concentration ratios / D b s by assuming that 
the ratio / O bs//o is equal to the ratio of the observed to the 
modeled values of 7 tot . By using the modeled weighted mean, 
lower limit, and upper limit of 7 tot we obtain the weighted 
mean, upper limit, and lower limit of / D b s . (Note that the lower 
limit of /tot gives the upper limit of / b s and vice versa.) Fi- 
nally, / b s is converted into the ratio of 5+ to total hydrogen 
number density (« 5 +/ w h) by multiplying by a factor of 1.1 
(assuming a helium to hydrogen number density ratio of 5%). 

Figure 13 shows the resulting ion concentration ratios as a 
function of height for the full range of model C data points. 
There is a hint of systematic radial increase at low heights. A 
similar radial increase would be predicted for ions that flow 
substantially faster than protons in the corona (by a factor of 
two) and then flow only ^10% faster than protons at 1 AU. 
Above 2 Rq, though, Figure 13 does not show any definitive 
radial trend. Taking account of the uncertainty limits, the data 
appear consistent with the 5+ ionization fraction being more 
or less "frozen in" (see, e.g., Hundhausen et al. 1968; Owocki 
et al. 1983; Ko et al. 1997). A linear least squares fit (using 
the logarithm of n 5+ / «h as the ordinate) is also shown, but 
the relatively high uncertainties at large heights preclude any 
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FIG. 13. — O s+ ion number density ratio (with respect to hydrogen) as a 
function of heliocentric distance for model C. Symbols show the weighted 
means of the reduced probability distributions, with styles the same as in 
Figure 1 . Vertical bars show the full range of parameter space with reduced 
probabilities greater than P[ CT . Also shown is a linear least-squares fit to the 
data points (dashed line) and empirical lower and upper limits as described 
in the text (gray region bounded by dotted lines). 

reliable interpretation of the slope. 

Figure 13 also shows a range of values that would have been 
expected from prior studies of both the oxygen abundance 
and the 5+ ionization fraction. The abundance ratio («o/"h) 
ranges from a relatively recent historical high of 8.5 x 10~ 4 
(Anders & Grevesse 1989) to the more recent — and some- 
what controversial — low of 4.6 x 10~ 4 (Asplund et al. 2004, 
2005; Grevesse et al. 2007). The ionization fraction (n s+ /no) 
was measured in situ by the SWICS instrument on Ulysses 
(Wimmer-Schweingruber et al. 1998) to be about 0.0058 in 
the fast solar wind. Models that include the freezing in of 
heavy ions have produced values for this ratio from 0.0035 
(Esser & Leer 1990) to about 0.005 (Chen et al. 2003). Al- 
though collisional ionization equilibrium is not expected to 
hold in the extended corona (for polar coronal holes), it is in- 
teresting to note that for T e = 10 6 K, both Arnaud & Rothen- 
flug (1985) and Mazzotta et al. (1998) give a ratio of about 
0.0045. This is similar to the above values, but it varies up and 
down by about a factor of 50% as T e is decreased or increased 
by only ±30%. We thus take tentative lower and upper limits 
for n 5+ 1 no of 0.003 and 0.006. The horizontal lines shown in 
Figure 1 3 were computed from the products of the two lower 
limits and the two upper limits given above for 

n^ = /no\/no 1 \ 

«H V"H/ V "O / 

The model C data points shown in Figure 1 3 have a mean 
value of n 5 +/ n H = 1-52 x 10~ 6 (taking a simple average) or 
1.39 x 10" 6 (taking the average of the logarithms). Perform- 
ing the same analysis using model R yielded values that were 
larger by about 5% (on average for all data points) to 20% 
(specifically for points at heights above ^3 Rq). The standard 
deviations of both sets of data points gave lower and upper 
bounds of approximately 8 x 10~ 7 and 2.4 x 10~ 6 that encom- 
pass the ±lcr range. The prior studies of oxygen abundance 
and O s+ ionization discussed above give somewhat higher val- 
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ues, which extend from 1.4 x 10~ 6 to 5.1 x 10~ 6 . Thus, both 
the model C and model R data points are in reasonably good 
agreement with the lower limit of the expected range, which 
gives some support for the recent low oxygen abundances of 
Asplund et al. (2004, 2005). 

5. DISCUSSION AND CONCLUSIONS 

The SOHO mission (Domingo et al. 1995) has made sig- 
nificant progress toward identifying and characterizing the 
processes that heat the corona and accelerate the solar wind 
(see also Fleck & Svestka 1997; Domingo 2002; Fleck 2004, 
2005). The results from the UVCS instrument regarding pref- 
erential heating and acceleration of heavy ions (i.e., 5+ ) have 
contributed in a major way to these advances in understanding 
over the past decade, and it is important to verify and confirm 
the key features of these results. Thus, this paper has ana- 
lyzed an expanded set of UVCS data from polar coronal holes 
at solar minimum with the goal of ascertaining whether ion 
temperature anisotropics are definitively present (as claimed 
by Kohl et al. 1997, 1998; Li et al. 1998; Cranmer et al. 1999; 
Antonucci et al. 2000; Zangrilli et al. 2002; Antonucci 2006; 
Telloni et al. 2007) or whether one can explain the obser- 
vations without such anisotropies (as claimed by Raouafi & 
Solanki 2004, 2006; Raouafi et al. 2007). These cases were 
exemplified by two sets of empirical models: one (model R) 
that was designed to replicate many of the conditions assumed 
by Raouafi & Solanki (2004, 2006), and one (model C) that 
used the same conditions as Cranmer et al. (1999). 

The main conclusion of this paper is that there remains 
strong evidence in favor of both preferential 5+ heating and 
acceleration and significant O s+ ion temperature anisotropy 
(in the sense 7]j_ > 7jii) above r w 2.1 Rq in coronal holes. 
More detailed conclusions, linked to the sections of the paper 
in which they were first discussed, are summarized as follows. 

1. It is important to search the full range of possible 
O s+ ion properties and not make arbitrary assumptions 
about, e.g., the ion outflow speed or the ion tempera- 
ture. It is clear from Figure 5b that if the comparison 
with observations is restricted to certain choices for the 
ion parameters, the resulting conclusions about the ion 
temperature anisotropy can be potentially misleading. 
(§4.1) 

2. The derived ion outflow speeds um and perpendicu- 
lar kinetic temperatures 7;j_ exhibit values similar to 
those reported by Kohl et al. (1998) and Cranmer et al. 
(1999), independent of the choices of electron density 
and flux-tube geometry. There is significant evidence 
for preferential ion heating and ion acceleration with 
respect to protons, although the radial rate of increase 
of Tj± may be slightly lower than that given by Cran- 
mer et al. (1999). The large values of Tj± appear to be 
due to true "thermal" motions and not unresolved wave 
motions. (§ 4.2) 

3. For heights above about 2.1 Rq, the models in this pa- 
per yielded higher probabilities of agreement with the 
UVCS observations for anisotropic velocity distribu- 
tions than for isotropic distributions. The UVCS obser- 
vations between the radii of 2.1 and 2.7 Rq were found 
to have probabilities of isotropy below about 10% (see 
Fig. 7), no matter what was assumed for the coronal 
electron density or flux-tube expansion (i.e., for either 



model R or model C). Even when using coronal proper- 
ties that seemed to maximize the probability of isotropy 
(e.g., model R), 78% of the UVCS data points exhibited 
probabilities of isotropy below our threshold one-sigma 
value of ~32%. (§ 4.3) 

4. The UVCS data at heights at and above 3 Rq can be 
used to put limits on the likelihood of strong 5+ tem- 
perature anisotropies. A key factor in discriminating 
between empirical models that either require or do not 
require a substantial anisotropy is the degree of exten- 
sion along the line of sight (LOS) of the emissivity. This 
extension is driven strongly by the rate of radial de- 
crease in the electron density. The relatively shallow 
slope of n e (r) used in model R (from eq. J8|) appears to 
be an "outlier" when compared to other observational 
and theoretical determinations of the electron density 
profile in coronal holes (see Fig. 3). Most other n e (r) 
curves exhibit a steeper radial decrease and thus a lesser 
degree of LOS extension for the O VI emissivities. Our 
model C, which utilized the empirical model parame- 
ters derived by Cranmer et al. (1999), had a represen- 
tative "steep" electron density profile and thus required 
a substantial 5+ temperature anisotropy to explain the 
UVCS observations above r w 3R Q . (§ 4.4) 

5. Models that exhibit enough LOS extension to reproduce 
the observed UVCS line profiles and intensities without 
a temperature anisotropy appear to require both (1) an 
electron density that decreases shallowly with increas- 
ing height, and (2) a highly superradial flux-tube geom- 
etry that projects a large fraction of the outflow speed 
vector into the LOS. Our model R, designed to be sim- 
ilar to the models used by Raouafi & Solanki (2004, 
2006), exhibited both of these conditions and thus had 
higher probabilities of an isotropic velocity distribution 
at heights above r w 3Rq. (§ 4.4) 

6. At the largest heights (r > 3Rq), the uncertainties in 
the existing UVCS measurements make difficult a firm 
determination of the anisotropy ratio. The analysis 
technique developed in this paper takes full account 
of these observational uncertainties. Future observa- 
tions with smaller observational uncertainties (see Fig. 
8) should yield correspondingly "sharper" probability 
distributions for the anisotropy ratio and thus better de- 
terminations of this quantity. (§ 4.3) 

7. Total intensities of the O VI AA1032, 1037 lines con- 
strain the ion concentration ratio n 5+ /«h to be approx- 
imately 1 .5 x 10~ 6 , with at least a factor of two range of 
uncertainty. If the freezing in of O s+ ions is considered 
to be relatively well understood, then this value implies 
a relatively low oxygen abundance in agreement with 
the recent downward revision of Asplund et al. (2004, 
2005). (§ 4.5) 

Because of existing observational uncertainties in the elec- 
tron density, flux tube geometry, and O VI line parameters 
such as Vi/e (the line width) and 1Z (the A 1032 to A 1037 inten- 
sity ratio), we cannot yet give "preferred" values for the 5+ 
anisotropy ratio 7jj_/7fii as a detailed function of height. Be- 
low r ps 2Rq, the observations are consistent with an isotropic 
velocity distribution. Between 2.2 and 2.7 Rq, the most prob- 
able anisotropy ratio appears to range between 2 and 10 (see 
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Fig. 6c). At heights around r w 3Rq the uncertainties are 
large, but there does seem to be evidence that the anisotropy 
ratio is likely to exceed 10 (see, e.g., Fig. 10b). The ratio must 
increase between 2 and 3 Rq, but we do not yet claim to know 
the exact rate of increase. 

New observations are required to make further progress. 
For example, as seen in Figure 3, there is still some disagree- 
ment about the radial dependence of electron density in polar 
coronal holes. Measurements of the white-light polarization 
brightness (pB) at solar minimum need to be made with lower 
uncertainties in the absolute radiometric calibration. Also, 
care must be taken to exclude time periods when high-latitude 
streamers may be contaminating the LOS in order to obtain 
a true mean electron density for a polar coronal hole. Ex- 
isting measurements of the superradial geometry (as traced 
by, e.g., polar plumes) tend to be limited by the fact that the 
shapes evident in LOS -integrated images are often assumed 
to be identical to the shapes of flux tubes in the plane of the 
sky. Better constraints on the flux-tube geometry could thus 
be made by using stereoscopy (Aschwanden 2005), tomogra- 
phy (e.g., Frazin et al. 2007), or other time-resolved rotational 
techniques (e.g., DeForest et al. 2001) to trace the full three- 
dimensional shapes of the plume-filled flux tubes. 

Improved ultraviolet spectroscopic measurements would 
greatly improve our ability to determine the plasma param- 
eters in coronal holes. We anticipate that the UVCS instru- 
ment will continue to observe polar coronal holes through the 
present solar minimum (2007-2008). We do not yet know 
whether the wide spread in line widths seen a decade ago 
(which exceeded the observational uncertainties) was due to 
a changing filling factor of polar plumes along the LOS or 
whether it may be connected to other kinds of time variabil- 
ity at the coronal base. An even wider range of variation in 
coronal hole properties was observed over the last solar cycle 
with UVCS (e.g., Miralles et al. 2006). These observations of 
how coronal holes evolve in size and latitude have helped to 
constrain the realm of possible parameter space of preferential 
ion heating and acceleration. 

There are also observations that cannot be made with 
UVCS that could greatly improve our understanding of ion 
energization in the solar wind acceleration region. For ex- 
ample, rather than having only 5+ (and, to a lesser extent, 



Mg 9+ ; see Kohl et al. 1999) in coronal holes, an instrument 
with greater sensitivity and a wider spectral range could sam- 
ple the velocity distributions of dozens of additional ions with 
a range of charges and masses. Obtaining the distribution of 
derived kinetic temperatures as a function of the ion charge- 
to-mass ratio Z/A would put a firm constraint on the shape of 
the power spectrum of cyclotron-resonant fluctuations (e.g., 
Hollweg 1999; Cranmer 2002b). A next-generation instru- 
ment with greater sensitivity may also be able to detect sub- 
tle departures from Gaussian line shapes that signal the pres- 
ence of specific non-Maxwellian distributions (e.g., Cranmer, 
1998,2001). 

The strong heating and acceleration of minor ions, as docu- 
mented by UWCS/SOHO, has provided significant insight into 
the physics of solar wind acceleration, but the basic chain of 
physical processes is still somewhat unclear. Many theoret- 
ical studies have attempted to trace this chain "backwards" 
from the known facts of kinetic ion energization to the prop- 
erties of, e.g., ion cyclotron resonant waves that can provide 
such energization naturally (see reviews by Hollweg & Isen- 
berg 2002; Marsch 2005; Kohl et al. 2006). Complementary 
progress has also been made in constraining the large-scale 
properties of the MHD fluctuations that may eventually cas- 
cade down to the microscopic kinetic scales (e.g., Verdini & 
Velli 2007; Cranmer et al. 2007). There are still many areas of 
disconnect, though, between our understanding of the macro- 
scopic MHD scales and the microscopic kinetic scales. Future 
theoretical work is expected to continue exploring how the 
combined state of turbulent fluctuations, wave-particle inter- 
actions, and species-dependent heating and acceleration can 
be produced and maintained. 
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FIG. 1 .— Collected UVCS polar coronal hole measurements of (a) O VI line widths V l/e , (b) ratio of O VI A1032 to O VI A1037 intensities, and (c) O VI 
A1032 line-integrated intensities, with symbols specifying the sources of the data (see labels for references). Error bars denote ±1ct observational uncertainties. 
Also shown (dotted lines) are the parameterized fits given by Cranmer et al. (1999). 



